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INTRODUCTION 

Measurement of vehicle dynamics during flight on test ranges has 
been approached using both ground-based and onboard instrumentation. 
Onboard platforms and accelerometers require telemetry systems or 
recorders resulting in relatively expensive instrumentation which must 
be considered expendable in most tests. Onboard solar aspect sensors 
have been developed for spinning missiles which can be used to infer 
yawing, pitching and rolling motion of flight vehicles. The most 
successful of these is the Yawsonde developed at the Ballistic Research ; 


Laboratories, ‘1*2) 


Many flight systems are tested using only ground- 
based instrumentation employing primarily photographic techniques. Data 
reduction in these cases is essentially all manual with attendant high 
cost and relatively long data reduction time. . 
During the past four years, two novel concepts for measurement of 
position and attitude have been investigated by the University of 
Wyoming in cooperation with the U.S. Army Missile Command. The first 
system investigated incorporated three ground-based laser transmitter/ 
detector stations and expendable retroreflecting elements located on the 
vehicles. This concept was investigated in considerable detail and the 
results of the studies are presented in reference 3. Research in the 
past year has been directed at a simplified system which requires only - 


two ground based stations. When compared to current ground-based 


photozraphic techniques for determining both position and attitude, 


Sa te 


XSite 


this system offers several advantages including: lower recurring costs; 
improved accuracy; automated data reduction; and applicability in low 
ambient light situations and in situations where the trajectory of the 


vehicle is not well established prior to flight. 


ne 
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SYSTEM CONCEPT 

The system includes two ground-based transmitter/detector tracking 
stations, each incorporating one pulsed and one continuous wave (CW) 
igsetice different frequencies. On board the vehicle, two different 
types of retroreflecting arrays are required. 

One type of array is composed of either conventional corner cubes, 
reflective tapes and or paints which have the property of retroreflecting 
a portion of the collimated incident beam back parallel to itself 
regardless of the orinetation of the reflective surface. This array 
forms a retroreflecting band located on the perimeter of the vehicle 
body at one axial position. Lilumination of and reflection by this 
retroreflecting band will then form the basis for a conventional laser 
radar tracking system. 

A series of roof type prisms as shown in Figure 1 form a second 
array which is also mounted on the vehicle's surface. This particular 
reflector array will be designated here the single plane corner reflector. 
A plane which passes through the center of all the reflecting surfaces 
of the roof type array and is also normal to all these surfaces will be 
called the retroreflection plane. Collimated light incident on the single 
plane corner reflector and contained in its retroreflection plane 
(B=90°) is reflected back to the source. Two of these single plane 
corner reflectors are mounted on the surface of the vehicle in such a 


manner that their retroreflection planes are skewed relative to each 


other and the roll axis of the vehicle. The mathematical analysis to 
follow will consider the special case where the retroreflection plane of 
one of the arrays contains the roll axis. An illustration of an 
instrumented vehicle is provided in Figure 2. 

4s the vehicle flies downrange, it is tracked with the laser radar 
for determining the vehicle position as a function of time while also 
positioning the CW laser to provide continuous CW illumination of the 
vehicle. During each revolution of the spinning vehicle, two CW laser 
pulses are returned to each of the two tracking stations. The time 
interval between the pulses returned to two separate tracking stations 
and between the two pulses returned to each station provide sufficient 
data for determination of the vehicle attitude. A mathematical description 


of this system is presented in the next section. 


MATHEMATICAL DESCRIPTION OF A SYSTEM 


In order to describe the system mathematically, two right-hand 
orthogonal coordinate systems are utilized: an earth-fixed cartesian 
system and a vehicle based cartesian system. The earth-fixed system 
is defined with the origin located at the launch site; positive Z-axis 
pointing in the vertical upward direction from the center of earth; 
positive Y-axis in the downrange direction; and positive X-axis in the 
crossrange direction. The vehicle based system is defined with the 
origin located at the vehicle center of gravity; w-axis coinciding with 
the vehicle roll axis with the positive direction pointing toward the 
nose; n-axis orientated perpendicular to w and parallel to the X-Y 
plane of the earth-fixed system; E~axis orientated perpendicular to 
n and w with its positive direction such that n, w, € form a right hand 
orthogonal system. 

The components of the position vectors of the ith ground station 
and the vehicle in the earth-fixed system are X 


i. Yo Z, and x Ya! 


Z respectively. Using the well-known transformation between the 


ny [cos(6_) sin(5,) 
wes sin(6,) ~sin(d,)cos(6,) cos($ 2) costs, ) 


cos (5, ) sin(é,)sin(é,) -sin(6é costs, ) 


6 


earth-fixed and vehicle-fixed systems one obtains the following rela- 
tionships for the coordinates of the ground stations in the vehicle- 


fixed coordinates: 
n, ™ cos (64) (X, - x) - sin(d,)(¥, - Y 
w, = sin(5,)cos(5,) (X, - x» + cos(6,)cos(6,) (XY, - Y 


+ sin(5,) (2, - Zz) (1) 


E, = ~sin(5,)sin(5,) (x, - x) - cos(5,)sin(5,)(¥, - Y. 


+ cos(6,)(Z, - Zz 


where 8) and 8, represent pitch and yaw respectively. Pitch is defined 
here to be the angle between the X-Y plane and the w axis whereas yaw 
is defined to be the angle between the YZ plane and the projection of 
the w axis into the X-Y plane as shown in Figure 3 and 4. 

(a) Time Relationship Between Return Pulses. 

To establish the time relationship between return pulses, we assume 
that the missile position, attitude and the roll rate do not change 
during the time interval between the pulse reception at the ground 
stations. Assume that at time ty there is a pulse returned to the 
first ground station from the skewed corner reflector as shown in Figures 
5 through Figure 7. At this moment if the position vectors of ground 
station 1 and 2 are Ny» Wy» by and Nos Wo» Eo » then it can be ob- 


served that, with a constant roll rate 2 for the flight vehicle, the 


time interval for reception of a second pulse from the second single 


plane array at the same station is 


1 
at, 7 ql2n - (a - 8.)] (2) 


The time interval between pulses returned from the first array to the 


two different ground stations is: 


1 , 
Ato, = 96) - &) (3) 


In terms of nw& coordinates, Eq. (2) can be written as 


w, tan(y) 
Mt), = 2n- a - gia Fa (4) 
nyt ey 
and Eq. (3) may be written as 
aT bo. 
Nbtos - aaa - paca Ey (5) 


Substituting Eq. (1) into Eqs. (4) and (5), it may be seen that the right 
sides of Eqs. (4) and (5) are functions of only two unknowns, that is, 
pitch (6) and yaw (8,). The relative positions are assumed to be 


known from tracking data. Since At and 2 can be determined 


11? 42) 
using data from the C.W. laser system, Eqs. (4) and (5) can in principle 
be solved for the pitch and yaw, using, for example, the Newton~ 


Raphson technique. 


(b) Mathematical Approach for Simulation 


The mathematical description presented in the previous section 


assumed that the vehicle was fixed in space, and that 2, 8) and 85 did 


not vary during the time intervals At,, and dt, In the situation 


8 


that arises during an actual flight test not only is the vehicle posi- 


tion changing but also 2, 5, and 6, may vary. One approach to data 


1 
reduction of flight test data would be to analyze the data as though 
the vehicle did meet the assumptions described in the previous section 
recognizing that some error could be introduced due to the vehicle 
dynamics. In order to gain some insight into the errors in the pitch 
and yaw angles determined using the static analysis a computer simu- 
lation a rocket trajectory was developed. The simulation provides the 
missile position and attitude as a function of time as well as the time 


at which pulses would be received from the two retroreflectors at the 


two ground stations. For this study a six degree of freedom trajectory 


(4) 


simulation originally developed by Harris 
10) 


» and later modified by 
Russel was again modified. The program developed by Russell was 
used to study the three ground station laser concept mentioned previously. 


Only those additional modifications required to simulate the two ground 


station - two single plane retroreflector concept are described here. 
In the original program, two right-hand orthogonal coordinate 


systems are defined: an earth-fixed coordinate system and a body~fixed 
coordinate system. The coordinate systems used in this thesis also 
include an earth-fixed and a body-fixed system, however, they are 
defined differently than the original systems. The relationships be- 


tween the coordinates originally defined and those used in this analysis 


may be seen in Figure 8. The body-fixed coordinates are obtained by 
rotating the earth-fixed system through three angles, A, uw, and v as 
shown in Figure 9. Other information which is needed in the mathe- 
matical approach of the simulation includes the direction cosines. 
These ~re defined in Figure 10. Then we have, from Figure 11 through 


Figure 13. 


8, = ~arcsin(£,) 


Lo 
8, ES arctan(>~) (6) 
1 


To determine the values of pitch (6,) and yaw (6,) at the moment when 
stations have return pulses, it is necessary to define an angle which 
is called the corner reflector angle (CRA). Suppose at time t < 0 

the skewed and straight corner reflectors are arranged on the rocket 
in the nw& coordinate system as shown in Figure 14. Imagine that the 
rocket is rotating in the clockwise direction. Then if at time = t', 
ground station 1 has a return pulse from the straight corner reflector, 
the rocket must be orientated in the nw— coordinates as shown in 
Figure 15. Considering this situation and referring to Figure 16 the 
corner reflector angle (CRA) can be defined as 


Ba " 
CRA = erctan(—) (7) 
5 


However, at this moment, the orientation for the position vector of 


the station 1 in nw coordinate system is 


I 
6, = arctan(—) 
1 ny 


10 


(8) 


This implies that the criteria to ensure pulses are returned from the 


straight corner reflector to station 1 is 


_ m3 1 
arctan(—) = arctan(—) 
ms, " 


This concept can also be extended to ground station 2. 


(9) 


Therefore, we 


can put the criteria to have pulses returned from the straight corner 


reflector in the general form: 


a, i 
arctan(—) = arctan(—) 
ny ny 


(10) 


To get a condition that station 1 has the signal returned from the 


skewed corner reflector, assume that at time = t, the rocket is 


eo! 


orientated in the n&w coordinate system as shown in Figure 15 and that 


station 1 has a pulse reception. Then, from the geometry, we have 


a- 6, = CRA + 6, 


Where 6, can be expressed as 


w,tany 


6, = arcsin[—37——,] 
(ny + €2)4 


From Eq. (7) through Eq. (12) we obtain 


(11) 


© (12) 


Re a eee ett 


+ 


11 


w, tany 
« - arcsin(—>—; 
(nj + a : 


my oy 
= arctan(—) + arctan(—) (13) 
3 a 


Eq. (13) is the criteria for pulse return at the station 1 from the 
skewed corner reflector. Similarly, for ith ground station to have 
pulses returned, the mathematical expression is the same as Eq. (13) 


except we replace 1 by i. 


(c) Sensitivity Analysis 

In the previous sections, when the equations for the returned pulses 
were set up to determine pitch and yaw, it was assumed that the missile 
position, attitude and roll rate did not change during the time inter- 
vale which were related to the actual reception of pulses at the two 
ground stations for the skewed and straight corner reflectors. As has 
been mentioned, the time intervals between returned pulses are influenced 
by the dynamic change of the missile position, attitude and roll rate. 
To determine the magnitude of the error which may be introduced by 
approximating the actual dynamic situation with a steady state model, 
we will define the following time intervals: bees is the observed time 
interval between the passing of the skewed and straight retroreflection 
planes through the ith ground station. Athy is the corresponding time 
interval predicted by the steady state model where the actual kinematic 


state occurring at the first received pulse is taken as the steady kinematic 


state of the missile. ate, ie the observed time interval between the 


CR eR at ee 


passing of the straight retroreflection plane through the ith and jth 
ground stations. Atty is the corresponding time interval predicted by 
the steady state model where the mean position, mean pitch and yaw, and 


mean rol] rate occurring between the two pulse times are used as the 


steady kinematic state of the missile, 


If the differences between the observed time intervals (at? 43° 


Att p and the corresponding time intervals predicted from the steady 


state model (act and att > are small, a linear variation between the 


ij 
two sets of time intervals in terms of kinematic differences may be 
appropriate. If this be the case, then the maximum error generated by 


the steady state model is approximately 


Ace, - ate cal (¢ - 4 y?# (14) ) 
if 44 = ao, t, "Kt, "Ke, 
aate 
O _ agh 14,2 - 2,35 
ats, ~ Aty, at My ye Cx,e, tow ] (15) 


where ty and to designate the time when the ith ground station has the 


actual pulse reception, represents the kth kinematic parameter 


"Et; 
evaluated at time ty» and ¢ represents fhe mean value_of the kinematic 
K,n 
dati, edt, 
parameter related to ates: iy ae and Be are ‘the 
- 
sensitivities of the time intervals atts bain ae. The peremsteh ty 


includes variables 6,, 5,, 2,, X_, Y. and Z. The various sensitivities 
1’ 2° kK OW Oe n 


The terms ( 


are classified into two categories for the skewed relfector and the 


straight reflector. 
(1) The Skewed System 


The various sensitivities for this system are pitch, yaw, altitude, 
range, crossrange, and roll rate. Mathematically they correspond to 


aati, dats, dAtys oat aAt aAt 


ii di ii 
a6, er az OY. OX and A respectively. Thus, if 
we use Eq. (4) and the parameter Pe to denote the variables 8) 85. 
\ 
; By. ee ear a oan . 


X,Y, and Z. We obtain 
en’ on ‘nm 


w an aft ow 


i i i i 
teva “aap? Si oes ae 
baud | RR ic eee Re RU ie (16) 
oP 2 2 2 2 4s 
K a 
a(n, + gy w,tan Y) 
and 
adt w,tany 
=a - 4 {a + arcsin[—}——> 7] - 20} (17) 
Q (ny + &5) 


(2) The Straight System 


The various sensitivities in this system have the same notation 


as in the skewed system except we replace At,, by At If PY indicates 


fi ij* 
the variables 8)» S55 x vy Zz then, from Eq. (5), the sensitivities 


are as follows: 


' 3g an de an 
ast "igpt - Sjapt "iggy - bi55 
ij 1 K K K K 
we, 7a th —3 7 ead Fee Gear cores Tae (18) 
K ni +e +E 
J j ames | i , 
and 
adt E gE 
ii. + farctan(—2) - arctan(—)] : (19) 
a2 2 ny nN 


The expressions for the sensitivities can be rewritten using equations 


(1) and the expressions presented below: 


rT ils (20) 


pa ~(X, > Xy)sin(8,)-(¥, - ¥ deos(6,) 


== ~ - X,)sin($,)cos(5,)-(¥, - ¥,,)cos(5,)cos(é,) 


- Zsin(s,) 


=" = -(X, - X,)cos($,)sin(5,)+(¥, - ¥,,)8in(6,)sin(6,) 


=y" = -cos(6 


pie -(X, > X,)ein(s,)sin(s,)-(Y, - ¥,)cos(6,)sin(5, ) 


+(Z, - Z,)cos(5,) 


- (X, - X,)cos(8)cos(4,)-(¥, - ¥.,)sin($,)cos(6,) 


2) 


— s sin($,) 


—s 0 


ag 


32 
a 


= sin(é,)sin(é,) 


4 
wr, = cos(§,)sin($,) 


« ~cos(6)) 


14 


(21) 


(22) 


(23) 


(24) 


© (28) 


(29) 
(30) 


(31) 
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—1. ~sin(6,)cos(6,) (32) 

vy." ~cos(S,)cos(é,) , (33) 
a 

du, 

oe -sin(6,) (34) 


(d) Error Analysis 


One approach to the data reduction would be to analyze the data as 
though the vehicle were fixed in space, with a fixed attitude (pitch and 
yaw) and a constant roll rate. The effect of vehicle dynamics would be 
to introduce an error into the values of pitch and yaw calculated from 
the measured time intervals. In order to obtain an estimate of the magnitude 
of the error that could ier daveodueea through this mechanism and also 
provide some insight with respect to which parameters play a dominant 
role in contributing to this error, the six degree of freedom simulation 
was exercised. The two ground stations were assumed to be located at 
3000, 0, O and 3000, 2000, 0. Parameters for the geometry of the skewed 
retroreflectors which had to be preassigned included a, the circumferential 
angle, and y, the inclination. Values of 60° and 15° were chosen for 


the respective quantities. 


The output from the simulation run included: 


a) For pulses retroreflected from the straight reflector 


1. The time of pulse reception at ground station 


number 1. (TIME1) 

The time of pulse reception at ground station 
number 2. (TIME2) 

The time difference between pulse reception at 

the two stations. (DTIME12) 

The time difference between DTIME12 and the no- 
minal time interval ati: (DTMF12) 

The missile position midway between the positions 
at the time of pulse receptions. (XM, YM, ZM) 

The average values of pitch, yaw and roll rate 
during the time interval DTIME12. (PITCH, YAW, ROLL). 
The difference in the average and the extreme 
values for missile position during the interval 
DTIME12. (DXM, DYM, DZM) 

The difference between the average and the extreme 
values of pitch, yaw and roll rate during the in- 
terval DTIME12. (DPITCH, DYAW, DROLL) 

The sensitivity of the time interval to: Pitch, 
yaw, crossrange, altitude and a time midway between 
TIME] and TIME2. (PITCH, YAWS, CROS, RANS, ALTS, 
ROLLS) . 

pulse retroreflected from the skewed reflector 

The ground station receiving the return pulse. (GSNO) 


The time at which a pulse is received at a particular 


station. (TIMES) 


16 


17 


3. Geometry related parameters described in the previous list. 
These results are presented in tabular form in the following 


section. 


pn ne re ee et en ee ne 
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RESULTS 


The six degree of freedom trajectory simulation was given the mass 
and aerodynamic properties of a tactical rocket (ARROW) recently tested 
at the U.S. army missile command. This vehicle was selected as repre- 
sentative of the type of vehicle for which a laser system might be used 
in a flight test program. Information which was derived from the simu- 
lation included: 

‘. Pulse reception time at the two ground stations from the 

two retroreflectors. 

2. Vehicle attitude, roll rate and position at the time of each 

pulse reception. 
This information was then used in a Newton-Raphson solution of Eqs. 
(4) and (5) to provide calculated values of pitch (5)) and yaw (55). It 
should be noted that average values of the geometric parameters were 
used in this solution which introduces the errors due to vehicle dynam- 
ics discussed previously. The pitch and yaw values obtained in this 
manner are presented in Table 1 where the values obtained from the simu- 
lation as well as the difference between the simulation and calculated 
results are also tabulated for comparison. 

The maximum difference between the simulation result and the cal- 
culated pitch is 0.81°, with an average difference of 0.15°, and a most 


probable difference of 0.03°. It is interesting to note that this 
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difference, or pitch error, is distributed along the downrange direc- 
tion with a maximum at the two cxtremes and a minimum near midrange. 
(see Figure 17) 

The maximum difference between the simulation result and the cal- 
culated yaw is 1.44° whereas the average difference is 0.18°. The most 
probable error is 0.04°. Referring to Figure 7 it may be seen that 
the yaw error is also range sensitive. In contrast to the pitch error 
distribution which shows a degradation in performance at each end of 
the range shown, the yaw error is large only at the extreme downrange 
positions. In order to investigate the sources of these errors the 
sensitivity analysis described in the previous chapter was initiated. 
The simulation was again used to generate the pulse reception times 
for the two reflectors and two ground stations. In addition the program 
was modified to calculate and print out the various sensitivities and 
perturbation terms. As described in the previous discussion the cal- 
culations were done in two groups that is for the straight retrore- 
flector with pulse reception at two different stations and for the 
skewed retroreflector system with pulse reception from the straight and 
skewed retroreflectors at the same station. The results are presented 
in Table 2 and 3 respectively. By multiplying the listed sensitivities 
by the appropriate perturbation terms the influence of each of the 
parameters on the time interval ati, or Ati may be obtained. The 
results of this operation are presented in abbreviated form in Table 4. 
In the interests of clarity only the order of magnitude of the parameters 
is indicated. It is immediately apparent that the error introduced by 


roll dynamics dominates the errors in At On the other hand roll, 


11° 


28 


position (range) and pitch dynamics all contribute to the errors in 
St), in a substantial way. Kelative magnitudes of the yaw and pitch 
contribution to the Aterrors may be thought of as yaw/pitch coupling 
in the At equation (Eqs. (4) and (5)). The Atio equation is weakly 
coupled to yaw throughout the first 1000 ft of flight and thereafter 
the magnitude of the pitch and yaw contributions become comparable. 

On the other hand the Ati, equation is weakly coupled to pitch beyond 
1000 feet and at ranges less than 1000 feet the pitch and yaw terms are 
comparable. A general observation which appears to hold for every 


parameter tabulated is the fact that the influence on the associated 


4t is strongest at the shortest ranges and longest ranges with a mini- 


mum in the region of 1000 feet. 
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EXPERIMENTAL PROGRAM 

(a)__Description of the Experiment 

The objective of the experimental program was to verify the system 
concept, in particular the mathematical formulation developed in the 
previous sections. It might be expected that problem areas not identified 
in the mathematical studies and earlier concept studies could be uncovered 
and investigated using the experimental apparatus. 

Two static Laser ground stations were fabricated, in addition to a 
model rocket equipped with two, 90° roof prisms. The ground stations 
were located on a line x = 30 ft. crossrange, and the rocket was systemat- 
ically positioned at various y and z coordinates along the line x = 0 in 
the earth fixed system. 

Each ground station incorporated a 3mw helium-neon laser operating 
at 6328 A,* a spatial filter and collimator. A mounting ring for the 
photodetectors was fixed to the collimator. Two holes were drilled 
through the ring and the photodetectors were inserted from the back and 
epoxied in place (see Figure 18). No auxiliary optics were used in the 
detection system since it was found that an adequate signal was obtained 
without this complication. A simple amplifier for the signal from the 
photodetectors was housed in a small aluminum chassis box which was 


epoxied to the top of the laser. 


*Metrogic Model 420 Laser 
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The mode] rocket consisted of a short, heavy-walled, alumigqum 
cylinder. At one axial position, two holes were bored through the wall 
of the cylinder. These holes were separated circumferentially by 
ninety degrees. One hole accomodated a mount for a straight reflector, 
and the other one accomodated a mount for a skewed reflector, allowing 
the angle y for the reflector to be adjusted to a particular value prior 
to conducting the experiment. The details of the design of the rocket 
may be seen in Figure 19. 

A steel drive shaft was attached to the back of the cylinder along 
the axis of the cylinder. The shaft was supported in cantilever fashion, 
with two ball bearings which were mounted in bearing blocks which were 
in turn attached to a base plate. An electric gear motor (9.61 rpm) 
was attached to the base plate and coupled to the drive shaft through 
a flexible coupling (see Figure 20). 

Both lasers were mounted on heavy duty metal tripods to facil ivace 
positioning and alignment. The model rocket was mounted on an L-bracket 
which was in turn attached to a rotary table. This arrangement provides 
a rather crude but rigid elevation over azimuth mount for the model. 

In order to measure the various time intervals, two digital timers 
were used.* The output signals of the photodetectors were amplified in 
the amplifier located on top of each laser. The outputs of these amplifiers 
were used as inputs to the timers to initiate and terminate counting. A 
schematic representation of the detection/timing circuit is shown in 


a re ee ee 


*Universal Timer Model 7370 
Hewlett Packard Timer Model 523DR 
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Fipure 2L.) The detection/timing etreuft was used to make three different 
types of measurements. Roll rate was measured using a signal from one 
laser ground station. The time interval for one complete revolution of 
the rocket was measured by starting and stopping the counter with consecutive 
er from a single photodetector. This measurement was checked with a 
nearly instantancous measurement of the roll rate which was obtained by 
atareing the counter with a pulse from a photodetector located in the 
mounting ring on top of the collimator, and stopping the counter with a 
pulse from a detector located in the same mounting ring directly below 
the center of the collimator. The time ‘interval between return pulses 
to the two laser scations from the straight reflector was measured by 
starting the counter with the pulse received at station 1, and stopping» 
the counter with the return pulse received at station 2. Simultaneously 
the time interval between pulse receptions from the straight and skewed 
reflectors at station 1 was measured. The same pulse which initiated 
counting in the previous measurement initiated counting on the second 
counter. This count was stopped however by a return pulse from the skewed 
reflector received at station 1. 

It is apparent that with only two counters available, all three 
intervals of interest could not be measured simultaneously. In view 
of the fact that the gear motor is synchronous, and therefore the roll 
rate would only change with the line frequency, the roll rate was 
measured separately and assumed to remain constant. Although this 
seems to be a reasonable assumption, it has been found that significant 


frequency fluctuations occur in the service provided by the University 


Power House. Therefore, the roll rate was measured before and after each 


attitude determination to check for constant roll rate. 


(>)__ Results of the Experimental Program 

Setting the angle y and measuring the various angles to a high 
degree of accuracy ultimately required the use of lasers and large optical 
levers. The entire system including ground stations, model, and retro- 
reflectors were eventually aligned to an angular accuracy of .05°. 
Experiments were conducted with the model located at two different positions 
downrange and pitch and yaw attitudes varying from ~20 to +20 degrees. 

The measured time intervals were deed to determine the pitch and 
yaw attitude using a Newton-Raphson technique as mentioned in the section 
entitled "Mathematical Description of the System'’. These tests indicate 
that, for the static system, the pitch and yaw can be determined to essentially 
the same accuracy as the experimental error associated with aligning 
the various components. Two demonstrations of this bench top system 


were conducted for army personnel. 


Pewee, - ones epee es eee 
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CONCLUS TONS 
In general, it may be concluded that both the analytical and 
experimental investigations of the skewed retroreflector two vround station 
Sy cesa made to date indicate that such a system could be a useful range 
instrumentation technique. The analytical studies indicate that the 
pitch error due to dynamics in the short range area arises from roll rate 


and range errors introduced into the At equation. In this region the 


12 
pitch and yaw solutions are weakly coupled and as a result the yaw 
solution remains quite good. Throughout: the midrange region the solution 
for both pitch and yaw appear to be accurate to within a tenth of a 
degree which could be adequate for range instrumentation. At the extreme 
downrange position, yaw and pitch solutions are more strongly coupled 
and large errors in both pitch and yaw are encountered which would be 
inadequate for a range instrumentation system. This phenomenon was, however, 
due to an error in the simulation which caused the rocket to experience 
extremely large pitch and yaw rates (1500°/sec) in the neighborhood 
1500 ft. downrange. These rates are unrealistically high. 

It shouldbe pointed out that while almost all aspects of the ARROW 
vehicle simulation are quite typical of tactical missiles and rockets 
which would be tested, the extremely large roll acceleration throughout 
the trajectory is atypical. Since the roll error is significant in both 
At equations, it is anticipated that the error terms would typically be 
smaller than those presented here. It is therefore concluded that the 
indicated degradation of system performance at 1500 ft. downrange ts 


not the result of an actual system limitation but rather the unrealistically 
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high angular rates generated in the simulation. 

Considerations such as ground station position optimization, retro- 
reflector geometry optimization, and alternate methods of data analysis 
have not been addressed. Further, work in these areas could be expected 
to produce a more accurate system. 

One limitation aot poluted out previously is the fact that the 
1500 feet range represents a maximum for retroreflection from the skewed 
reflector to ground station number 1. It is now apparent that somewhere 
near midrange an equation for Atos should have been substituted for Atay 


which could have allowed data reduction out to 2000 feet and could 


possibly result in imporved pitch and yaw solutions between 1300 feet 


and 2000 feet downrange. It should be noted that the yaw accuracy obtained 


is better than that associated with a similar system employing three 
ground stations and one single plane corner reflector onboard the 
ventete 

In summary, the two station concept offers all of the advantages 
orginally anticipated for the three station system with the additional 
economic advantage of eliminating one ground station. It appears that 
yaw accuracy will be better than that associated with the three station 
system. Additional study covering ground station position optimization, 
retroreflector geometry optimization, and data analysis solutions are 


all required prior to implementing the concept as a range instrumentation 


technique. 
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center line of skewed 
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straight corner 
reflector 


VEHICLE END VIEW 


FIGURE 5 - CORNER REFLECTOR GEOMETRY RELATIVE TO VEHICLE 
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FIGURE 6 ~ GEOMETRY FOR ANALYSIS 
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FIGURE 7 - GEOMETRY FOR ANALYSIS 
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(a) EARTH-FIXED 
COORDINATE SYSTEM 


571 FT 


XX 
(b) EARTH-FIXED COORDINATE SYSTEM 
AS DEFINED IN THE SIMULATION 
PROGRAM 
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FIGURE 8 - RELATIONSHIPS BETWEEN TWO EARTH-FIXED 


COORDINATE SYSTEMS 
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FIGURE 9 ~ RELATIONSHIP BETWEEN THE BODY FIXED 


AND THE EARTH-FIXED COORDINATE SYSTEM 


A rotation about the ZZ-axis through the angle \ produces the 
X'Y'Z' -coordinate system where Z' = ZZ. Then a rotation about 
Y’ through the angle u produces the X"Y"Z" ~-coordinate system 
where Y" = Y'. Finally, a rotation about the X" -axis through 
the angle v produces the XnY 7B -coordinate system. 
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FIGURE 10 - DEFINITION OF 


DIRECTION COSINES 
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FIGURE 11 - DIFFERENCE BETWEEN nw AND THE 


BODY~FIXED COORDINATE SYSTEM 


n-axis always parallel to X-Y plane, Yg and Zg axis rotate at 
the same angular velocity as the vehicle does. 
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FIGURE 12 - PITCH ANGLE IN TERMS OF 
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XX = YY 


YY = X 


FIGURE 13 - YAW ANGLE IN TERMS OF DIRECTION COSINES 
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skewed corner reflector 


straight corner reflector 


FIGURE 14 - GEOMETRY OF CORNER REFLECTOR AT t = 0 
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FIGURE 15 - GEOMETRY RELATING TO THE CORNER 


REFLECTOR ANGLE 
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FIGURE 16 - CORNER REFLECTOR ANGLE IN TERMS 


OF DIRECTION COSINES 
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FIGURE 17 - PITCH AND YAW OF ERROR DISTRIBUTION 


WITH DOWNRANGE DISTANCE 
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